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SUMMARY 


A general nomographic solution for the induced velocities and 
wake skew angle is presented for use with VTOL-STOL systems. 


INTRODUCTION 


VTOL-STOL aircraft are characterized in general by the fact that 
in some portion of their flight envelope the wake is sharply inclined 
to the free stream. Under such conditions, the usual small- angle 
assumptions used in determining the induced velocities, and conse- 
quently, the power required, are no longer valid. Indeed, the use of 
small- angle assumptions leads to such anomalous results as infinite 
induced velocities and required power in the extreme case of hovering. 

The aforementioned difficulties may be avoided by a more complete 
examination of the horizontal and vertical momentum imparted to the air 
by the aircraft at low speeds. The resulting equation is a quartic in 
the induced velocity, and, as such, is difficult to apply. On the other 
hand, this quartic can be solved in its most general terms and the 
resulting solution then can be derived and presented in the form of a 
chart, or nomograph, from which the required induced velocities may be 
read directly. This paper presents such a chart. 

The problem of estimating the induced velocities becomes of addi- 
tional importance since the induced velocities combine with the forward 
velocity to determine the wake skew angle upon which, for example, wind- 
tunnel corrections and ground effect have been shown to depend. (See 
refs. 1 and 2.) The wake skew angle also can be read directly from the 
nomograph presented herein. 



SYMBOLS 


momentum area of aerodynamic force- generating system, sq ft 
induced drag, or horizontal component of force, lb 
lift, lb 

ratio of initial to final induced velocity 
free- stream dynamic pressure, ipV^, lb/ sq ft 

slipstream dynamic pressure, ^p^nuQ + V)^, lb/sq ft 
thrust of a rotor, lb 

horizontal induced velocity, at the force- generating system, 
required to produce a given horizontal force, positive rear- 
ward, ft /sec 

free- stream velocity, ft/sec 

resultant velocity at force- generating system, ft/sec 

vertical induced velocity, at the force- generating system, 
required to produce a given vertical force, positive upward, 
ft/sec 

hovering induced velocity, vertical induced velocity when 
V = 0 and % = 0 ( see eqs. ( 8 ) and ( 23 )), positive 

upward, ft /sec 

lift and drag axes, perpendicular and parallel, respectively, 
to the horizontal or free-stream direction (see fig. l) 

angle of attack of tip-path plane, positive when leading edge 
of disk is high, deg 

net wake deflection angle, angle between the wake and the 

X-axis, or free-stream direction, 90° - X, deg (see fig. l) 

mean induced velocity, normal to the tip-path plane of a rotor 
positive downward, ft/sec 
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hovering mean- induced velocity, normal to the tip-path plane 
1 of a rotor, positive downward, ft/sec 

p mass density of air, slugs /cu ft 

X wake skew angle, angle between the wake and the negative 

Z-axis, positive when measured rearward, deg (see fig. l) 


THEORY 


Induced Velocities 

Because of its lift and drag forces, the machine creates the veloc- 
ities uq and wq at its own location and the velocities nu^ and uwq 

in the far wake. The force and velocity vectors are shown in figure 1 . 
From this figure, it may be seen that the resultant velocity Vr at the 

force- generating device is 


Vr = ]f(v + uq ) 2 + (-wq ) 2 (l) 

Furthermore, the lift and induced drag are 

L = PAn v R( _nw o) (2) 

D i = PV^-nuo) (3) 

where ^ is the area of the cylinder of air affected by a wing, the 

actual area of a propeller or rotor, or the duct- exit area of a ducted 
fan. 


Notice that a forward thrust is merely considered equivalent to a 
negative drag in equation ( 3 ). 


Dividing equation ( 3 ) by equation (2) yields 

Di = uo 

L w 0 


(M 


provided, of course, that L 


is not zero. 



Now substitute equation (4) into equation (l) to obtain 



Now, if the force- generating device could hover with the same lift, 
the same momentum area, zero forward speed, and zero induced drag, the 
vertical induced velocity in hovering w h would be found from equa- 
tions (6) and ( 7 ) as 


w h = ± 



(8) 


where the negative sign corresponds to a positive or upward lift and 
the positive sign corresponds to a negative or downward lift. Note 
that w^ is purely fictitious for a wing, but that it can still be 

calculated for use as a reference velocity. 

Combining equations (6), (j), and (8) yields 

( M 2 = 2^0 
V w h/ v 



( 9 ) 


or, finally, squaring both sides of equation ( 9 ) 
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Equation (10) is essentially the same as the equivalent expression 
derived for a rotor in reference 3 except that in the present case equa- 
tion (10) gives the vertical component of induced velocity in terms of 
an arbitrary value of the induced drag- lift ratio, whereas reference 3 
solves for a total induced velocity which is assumed to be normal to 
the plane of the rotor. 

Equation (10) is, of course, a quartic in Wq and, as such, is 

difficult to solve for any one given case. Points for use in the con- 
struction of a general nomographic solution, however, may be obtained 
rapidly by substituting a series of values for V/wq into equation (10), 

solving for wq/ w^, and then noting that 


V_ = JL 

W h W 0 W h 


( 11 ) 


The results of this calculation can then be plotted with w q/w^ 
as a function of V/ w^ for specified values of D^/L. Such a chart is 

presented in figure 2; however, a full discussion of it will be deferred 
to a later section of this paper. 

With wq obtained from a general nomograph of equations (10) 
and (11), uq can be obtained from equation (4). 


Wake Skew Angle 

The remaining item of interest is the wake skew angle, or its com- 
plement, the angle by which the wake is deflected from the horizontal 
or free- stream direction. This angle may be obtained by examining fig- 
ure 1 which indicates that 


V + UQ 
-Wq 


tan X 


V Dj 
wq ” T 


( 12 ) 


so that equation (9) may be rewritten as 


w, 


0 


or, solving equation (13) for X 


(13) 


lxl.c-1 S' 


w. 


(14) 


where , from equation (12), X is positive if — — > — and negative 


if 4-<3t 


-W 0 


-Vr 


The net wake deflection angle, being the complement of X, is then 
immediately given as 



(15) 


and, similarly, 0 n is in the first quadrant if 

V Di 

second quadrant if < . 

-w 0 L 



and in the 


The Case L = 0 

When L is zero, the foregoing treatment is not necessary. Note 
that for L = 0, wq = 0, so that equation (l) becomes 

V R = V + UQ (16) 

Equation (3) may then be solved immediately in the form 



It is of interest to square both sides of equation (17) and then 
multiply both sides by 2p to obtain 
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ip(2uo + v) 2 =ipy2.§^ ( 18 ) 

which for the large class of devices characterized by n = 2 (i.e., 

wings, rotors, unshrouded propellers) becomes 

<1s = * - U9) 

Equation (19) is, of course, quite familiar since it yields the 
slipstream dynamic pressure used in nondimensionalizing many VTOL test 
results. (Notice that the thrust of a propeller at zero angle of attack 
is equal to -D^.) It is evident that equation (19) holds for all wings, 
rotors, and propellers where n is equal to 2. On the other hand, con- 
sider, for example, a ducted fan where n is not equal to 2. In this 
case, equation (17) cannot be simplified to the form of equation (19). 
Note that, in general, 


Is - |p( n,J 0 + v) 2 

/ 9 - ( 20 ) 

It is, of course, obvious for the case of L = 0 that X = 90° 
and 0 n = 0°. (See fig. 1.) 

Equivalence of Reference 3 

Reference 3 derives an induced- velocity equation which corresponds 
to equation (10) for the special case of a rotor. Since the work of 
reference 3 has been confirmed experimentally, both with regard to 
induced velocity (ref. 4) and wake skew angle (ref. 5) it is advisable 
to examine the equi valence of reference 3 and the present analysis. 

Reference 3 assumes a thrust force T on the rotor which is nor- 
mal to the rotor disk and a corresponding induced velocity v also 
normal to the disk. Thus, from figure l(b) 


— = tan a 
L 


( 21 ) 
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v 0 



2 

+ U(f 


= w, 


°i 


i + 



-vp 

cos a 


( 22 ) 


v h = 


- v h 
cos a 


(23) 


Substituting the values of w q, w^, and %/L as obtained from 

equations (21), (22), and (23) into equation (10), expanding, and then 
simplifying leads to 



(24) 


which, noting that a in this paper is identical to -a in reference 3, 
is identical to the expression derived in reference 3- 

In view of the complete equivalence of reference 3 and the present 
analysis, the experimental confirmation of references k and 5 may be 
considered to apply to the present analysis as well. 


PRESENTATION OF RESULTS 


A general nomograph for solution of the momentum equation, based 
on the preceding equations, is presented in figure 2, This figure is 
primarily a set of curves giving wq/w^ as a function of Vy w^ for 

specified values of the induced drag-lift ratio. As such Wq/w^ may 
be read directly on the ordinate for any given value of V/w^. It is 
understood, of course, that the velocity w^ is obtained from 
equation (8). 

In consequence of equation (ll), lines of constant V/wq are 

radial lines in figure 2. A number of such lines are shown. These 
may be used as a guide in determining the values of V/wq without 
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recourse to equation (ll). In this regard, for V/wq ^ 6, a line may 

be drawn through both the origin and any point determined by specified 
values of V/w^ and Dp/L. The required value of V/wq may then be 

read directly from the V/w^ scale at the intersection of the radial 
line and the top of the figure (wq/w^ = l.o). 

According to equations (l4) and (l5), X and 0 n are functions 
of wq/w^ only. Therefore scales have been added to figure 2 in order 

to allow these quantities to be read directly. These scales are 
extremely compressed for wake angles near the vertical. For greater 
ease of reading, figure 3 has been prepared. This figure presents 
directly the wake angles as a function of wq/w^. Note in particular 

the signs and quadrants of these angles as pointed out following equa- 
tions (l4) and (15). 

The horizontal induced velocity uq can be obtained by using 
w 0 /w h , as obtained from figure 2, in equation (4). In the case L = 0, 
uq may be obtained from equation (17). 

The same chart may be used in rearward flight as well provided that 
V/ W h, V/Wq, X, and Dp/ L are all simultaneously multiplied by -1. 


DISCUSSION 


It will be observed in figure 2 that, for all cases of forward 
thrust (Di/L < 0), the vertical induced velocity Wq decreases con- 
tinuously as the forward velocity increases. Note that when Dp/L < 0 

there is a forward- directed component of force and, consequently, a 
rearward- directed horizontal induced velocity even at zero forward 
speed. The superposition of a forward velocity always adds to this 
horizontal velocity, increasing the mass flow through the system and 
reducing the required vertical induced velocity for a given lift. 

On the other hand, when there is a drag or negative thrust 
(Dp/L >0), the vertical induced velocity increases until wq/w^ is 

equal to one and then decreases thereafter. The reason for this 
behavior will become evident upon consideration of the mass flow 
through the system. In this case, when Dp/L > 0, the horizontal 
force is rearward and the net horizontal velocity is the difference 
between the forward velocity and the horizontal induced velocity. 

Thus, as the forward speed is increased from zero, the mass flow is 
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decreased and the vertical induced velocity required for a given lift 
is increased. At a forward speed of V/w^ = V/wq - -D^/L, the forward 

velocity exactly cancels the horizontal component of induced velocity. 
Thus, at this speed, the mass flow has its minimum value and the 
required vertical induced velocity has its maximum value. Further 
increases in forward speed increase the mass flow so that the vertical 
induced velocity decreases thereafter. 

At very large values of D^/L, the curves of wo/w h against V/w h 
become triple valued in the range of values near V/w^ of -Dj_/L. For 

these conditions, it is likely that no true wake, as postulated by 
momentum theory, exists. This situation is similar to the vortex ring 
state of a propeller or rotor (ref. 6) and probably also to the sepa- 
rated flow of a flat plate at high angles of attack. 

It should be noted in particular that the induced velocities and 
the wake skew angle of the present analysis are always defined in rela- 
tion to the lift and drag axes of the aircraft. This is not always the 
case when applying the results of the analysis and suitable adjustments 
must be made in certain cases. For example, reference 7 gives the flow 
field of rotors or propellers as a function of wake skew angle where 
the skew angle is defined in relation to the tip- path- plane axes. In 
this case, X of reference 7 is equal to X + a from the present 
paper. 


CONCLUDING REMARKS 


A nomographic solution for the induced velocities and wake skew 
angle of an arbitrary aerodynamic force- generating system has been given. 
These quantities may be obtained directly from a chart which requires a 
knowledge of only the lift and induced drag of the system. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., February l4, 1961 . 
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